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I. DUAL-PHASE CYCLE 


The dual-phase nozzle is a key element in the dual-phase 


engine concept. This nozzle will be studied in detail, but 


first a review of the dual-phase cycle will be carried out | 
from information obtained from References 1 and 2. 
Reference 1 describes the dual-phase engine. It is a 
new concept which operates on a mixture of two fluids cr jo 
| phases of one fluid. This cycle employs a high-torque/i: - 
rpm impulse turbine which eliminates the requirement for 
speed reduction needed in a conventional steam turbine. 
This allows for direct drive of a ship's propeller with no 
gear reduction. This is ideally suited for marine propul- 
sion, reducing weight, noise, and volume. 
The dual-phase concept is a modified Rankine cycle with 


subtle differences of extreme importance. In a normal steam 


turbine the working fluid enters a turbine through a nezzle 


where the kinetic energy of the working fluic is converted 


to a mechanical form. The dual-phase system introduces a 
second fluid prior to entry into the nozzie. This fluid is of 
higher vapor pressure than the steam and therefore remains 
7 in the liquid state throughout the cycle. Section A will ., 
describe this two-component cycle while section 2 will do 


the same for a single-component system. The dual-phase 


system can be divided into the two groups illustrated in 
Figure 1. Single-component flow can be further divided into 


three categories. 


A. TWO-COMPONEN 
A two-component mixture is one in which the lew vapor 
pressure liquid and a high vapor pressure liguid are of dif- 


ferent chemical compounds. Some fluid combinations which 


have been considered are steam-krytox, steam-caloria, steam- 
lead, bismuth eutectic, and dow-therminol. A schematic flow 
diagram and process representation on the T-S diagram are 


shown in Figures 2 and 3 for the two-phase engine cycle using 


> 


a "two-component" mixture. The liquid phase is lith 


bP 


um car- 
Donate and vapor phase is steam. To illustrate the overall 
advantages of the two-phase engine cycle a discussion on the 
theory of operation will be presented using a two-compcnent 
Mixture and Figures 2 and 3. 

The major component of the dual-phase system is the nozzle. 
A mixer area is located prior to the nozzle inlet. A high 
vapor pressure liquid is placed in contact with the low vapor 


pressure liquid in this area. A nigh pressure vapor liquid 


greater than temperature of the water, heat is transferred 
to the water causing it to vaporize to point 1. Figure 4 
illustrates the temperature and state point from the inlet 


af the mixture area *o the nezzle exit. This mixture is 
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Figure 3. Dual-~Phase Two-Component T-S Diagram 


mmm 


OFJeuMsYyoS suTquny estnduy pynbry *h asanBTy 


(2/A AMDO13A) ae 
"dO0D5 ainol 


(0¥3Z ALIDO13A 31M1OSAV 
“WA ALMDOTIA 3AILV139) 
“INO ainNDN! 


INIGINL GINO 


a) 
ee 
= 
Sa 
an 


JAVHS INdLNO Y3IMOd 


aN 


(A ALIDO13A WNYQ) 
“YOLVYVdIS AYVLOY 
ONINNIdS-3383 


hess eamdnumtenteeny ee tee teehee tee) 


O) 
4 Uj 


(bh rcbntrbhheaheahethtase 


(A ALIDOTSA 13f) 
‘7IZZON 3SVHd-OML 


14 


eae 
nn eee ee ee 


expanded to a low pressure at the nozzle exit. The expand- 
ing vapor transfers momentum to the liquid droplets, while 
the droplets transfer heat to the vapor. This tends to ap- 
aproach isothermal expansion through the nozzle, point 2. 
After leaving the nozzle, the liquid droplets are separated 
from the vapor by a rotary separator located on the turbine. 
The rotating drum turns at close to the nozzle's exit vel- 
ocity. Hence there is very little friction loss between the 
separated liquid and the drum wall. The centrifugal accel- 
eration produces a very clean separation of vapor and liquid. 
The liquid traveling with the rim of the rotating drum trans- 
fers axially through holes in the separator disc and subse- 
quently enters the liquid turbine proper. The kinetic 
energy of the liquid is converted to shaft horsepower by the 
liquid turbine. Figure 5 illustrates the dual-phase liquid 
impulse turbine assembly. The vapor remains superheated as 
a consequence of the heat transfer from the liquid. The 
vapor flows from the separator, point 3, through a regenera- 
tor where heat is transferred to the condensate. The vapor 
is condensed, point $; pumped to nozzle pressure, point 6; 
and passed through the regenerator for heating. Heat is 
added Detween point 6 and point 1 in twos methods. The re- 
generator adds heat to the condensate by using the steam 
from the turbine rotor, point 7. The remainder of heat,‘ 
point 7 to 1, is added by the heated liquid mixed with she 
condensate in the mixer. The water is vaporized by direct- 


contact heat transfer. 


where 1, = Temp of liquid entering mixer 


T, = Temp of water entering mixer 


T, = Temp of steam leaving the nozzle 
T = Temp of liquid leaving the nozzle 
L 

2 


Figure §. Temperature & State Point Diagram 


for the Mixture & Nozzle 


In the nozzle most of the thermal energy of the steam 
is converted to kinetic energy of the liquid droplets. 
This acceleration of the liquid by the vapor in the two- 
phase nozzle provides the kinetic energy to drive the 
liquid impulse turbine. The liquid velocities involved 
are relatively low as compared to velocity of the vapor. 
Thus, the output of the impulse hydraulic turbine will be 
high torque/low rpm. This conversion of the liquid kinetic 
energy to shaft power at high torque with low rpm appears 


to have direct application to naval propulsion. 


B. ONE-COMPONENT 

A one-component system is one in which the working fluid 
is of the same chemical compound. One of the simplest dual- 
phase one-component systems is illustrated in Figures 6 and 
7. The working fluid is heated to saturation temperature by 
some type of heat source. Heat sources applicable to this 
case are geothermal power olants, engine exhaust, industrial 
Wast~-heat recovery, and bottcming cycles for steam and gas 
turbine plants. This working fluid, at saturated liquid con- 
ditions, with small amounts of vapor is expanded through a 
two-phase nozzle. As the expansion process takes place ,the 
liquid partially vaporizes and accelerates the remaining 
liquid phase in the nozzle. The dual-phase mixture enters 
the rotary 3eparator and the same process occurs as mentioned 


in section A. Since the liquid phase is of a much higher 
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Figure 7. Dual-Phase Single-Ccmronent T-3 Diagram 
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density than the vapor, the liquid velocities are relatively 
low as compared to the velocity of the vapor. Thus the out- 
put of the turbine will again have high torque at low rpm. 
After leaving the turbine the vapor mixture is condensed and 
the condensate is pumped back to the heat source. The cycle 
is shown on a T-S diagram in Figure 6. The state points are 
numbered to correspond with Figure 5. The dual-phase nozzle 
expansion takes the fluid from a saturated liquid, point 1, 
to a dual-phase flow, point 2. The flow is deceierated in 
the rotor; condensed, point 4; and pumped back to nozzle 
inlet pressure at point 5. The liquid is then reheated by 
the source fluid to point l. 

Another application of the one-component two-phase cycle 
is the wet-to-dry cycle. If the initial temperature cof the 
working fluid is sufficiently high and the saturation curve 
nas a positive saturated liquid slope the working fluid can 


rm 
£ 


be exvpanded to dry vapor. igure 3 is the T- diagram for 
a wet-to-dry cycle. The fluid is expanded from saturated 
liquid at point 1 to saturated vapor at point 2... The vapor 
drives animpulse rotor and leaves the rotor slightly super- 


heated at point 3. The vapor is condensed to point 4 and 


pumped back to the nozzle inlet pressure at point 5. 


C. ADYANTAGES 
The advantage of the dual-phase cycle with respzect to 


Marine application is the ability to achieve low shaft speed 
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igure 8. Wet-to-Dry T-S Diagram 


in a small turbine engine. For example a steam turbine of 
150-kw shaft horsepower using the temperature given in 
Figure 2 would have a speed of about 60,000 rpm. A compar- 
able two-phase turbine would operate at approximately 10,000 
rpm. There is also an efficiency advantage. At temperature 
corresponding to Figure 2, a steam Rankine cycle would have 
an efficiency of approximately 23% where-as the two-phase 
cycle would have an efficiency of 37%. This is assuming equal 
turbine efficiencies. The two-phase cycle also allows for 
control of turbine speed because the vapor/liquid mixture 
ratio can be varied to change the nozzle exit velocity. This 
is a capability unavailable in a conventional steam turbine. 
Both of the dual-phase concepts can be thought of as a 
form of a regenerated Rankine cycle. The dual-phase cycle 
Dy control of liauid/vapor mixture ratio enhances the overall 
power system controllability. The T-S relationship for a dual- 
phase two-component engine cycle compared to a Rankine cycle 


—™s 
. 


is shown in Figure 9. 
Two design studies References Land 2, have shown potential 
advantages in the two-phase engine cycle as compared to the 


conventional Rankine cycle fcr marine propulsion. The fol- 


lowing advantages were noted: 


r 


1. High Efficiency - Full lead output power performance 
gains ranging from 20 to 59 percent was found. 
2. Direct Drive - Direct drive at speeds ranging from 


90-4500 rpm was found possible with a single-stage turbine. 


dual-phase 
expansion 


f 


2,3 


ideal steam 
cycle expansion 


TEMPERATURE 


ENTROPY 


Figure 3, T-S5 Comparison of Dual-Phase Two-Component 
Cycle & Rankine Cycle 
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3. Reduced Volume - Volume reduction of 30 percent 


were estimated. 


4u, High Part-Load Efficiency - Variable mass ratio 
enabled part-load (cruise) efficiency gain of as much as 


1090 percent. 
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IIT. DUAL-PHASE NOZZLE THEORY 


The flow phenomenon of a two-phase mixture has been ana- 
lvzed in Reference 1. It is repeated as follows. The problem 
is illustrated in Figure 10. A spatially uniform two-component 
mixture of liquid drops and gas enters a nozzle at high pres- 
sure and low velocity and expands to low pressure and high 
velocity. The objective of the analysis is to determine, for 


a specified pressure the drop diameter D and the tempera- 


tures T_ and Tho velocities VY, and Vy and flow rates ms and 
—) 


m, ef the gas and liquid phases, respectively, at each station 
~s 
in the nozzle given the initial values of D, T_, T,, V_,%,, 
sS ~ 5 ~ 


the total flow rate, and the properties of the fluids. 


Figure 10. Sual-Phase Nezzie Flow 


Seometry & Nomenclature 


The five relations employed to compute the five unknowns 


D, Ty, Tyr Vos 


mixture, (2) the energy equation for the mixture, (3) the 


aM and V;, are (1) the momentum equation for the 
drop drag equation, (4) the drop heat transfer equation, and 
(§) the drop breakup criterion. Solubility and vapor pres- 
sure relations provide the flow rate ratio mo / my 
A. ASSUMPTIONS 

The assumptions employed in the two-component analysis 
are as follows: 

1. The liquid is uniformly dispersed as spherical drops 

all of the same diameter. 


2. The drops break up to limit the Weber number to 6. 


Ww 


. There are no external forces acting on the two-phase 

mixture other than pressure and wall shear, and 

there is no heat transfer to or from the mixture. : 
4, The flow is one-dimensional. 
5. The drops are large enough for the surface curvature 


to Nave negligible effect on the vapor pressure o 


the liguid and for the surface energy to be negligible. 


om 
rj 
on 
o 


drops are isothermal. 
?, The gas mixture obeys the additive-pressure law. 


3. The partial pressure of the predominantly liquid 


# 
[ey 


component in the liquid is given by Raoult's Law. 


9, The concentration of the predominantly gaseous com- 


ponent in the iiauiz is given py Yenry's Law. 


pre ET a ee ee. 


10. The volume of the liquid solution is equal to the 

sum of the volumes of the pure licuids. 

Assumption 1 restricts the analysis to nozzles having 
Spatially uniform injection of the liquid into the gas and 
Operating at gas-to-liquid volume ratios greater than unity. 
Assumption 2, the drop breakup criterion, states that drop 

“pe 


diameter is limited to a value D for which Se = ae w/e 
— 


= 6. Thus 


where O5 is the gas density, We is the slip velocity V, - Vo 


and is the liquid surface tension. The form of =a. (1) is 
physically reasonable in that the Weber number is propor- 
tional to the ratio of stagnation pressure ee 5/2 to sur- 
face tension pressure 4V/D. Hence, a drop would be expected 
to flatten and breakup at a sufficiently high value of “2. 
This has been verified experimentally and the critical Weber 
number found to be 3, within a factor of about two. sal 
additional restriction is that for actual Sreakup to occur, 
the time spent at a Weber number exceeding 5 must be longer 
than the natural period of oscillation of the drop, 


mp, D°/ prety 


4, where Oy is the density of the liquid. 

This requirement is met only in two-phase nozzles longer 
than about 10 in. and Assumption 2 may cause the analysis to 
overestimate the exit velocity by increasing amounts as the 


nozcle length decreases below 19 in. 
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Assumption 3 excludes magnetohydrodynamic and mechanical 


body forces. The exclusion of wall heat transfer is correct 
for the insulated nozzles of interest for power systems. In 
addition, the relatively high velocity results in short 
residence lines in the nozzle proper. 

Assumption 4 1s closely met in practical nozzles since 
good performance requires small wall angles, large throat 
radius of curvature, and uniformly distributed injection of 
the fluids at the nozzle entrance. 

Assumption $ is valid for the drop sizes of 9.9C1 to 
0.910 in. produced by the Eq. (1) breakup criterion. Assump- 
tion 6 is valid because of the rapid internal circulation 


in drops. Assumption 7 introduces negligible error in most 


cases of practical interest since the vatror pressure of the 
liguid is small and needs only to be evaluated approximately. 
A.tumptions 3, 9, and 10 are either valid, or cause little 


error, for fluids of low miscibility, which are the fluids of 


interest. 


3. DERIVATION OF EQUATIONS FOR FREE-STREAM FLOW 
1. Contiauity 
Referring to rigure 1, the nozzle flow area A is 


equal to the gas flow area moO ov, plus the liquid flow 
ea tt \ 
area mh, /0,%,- Thus . 
> 1 Ld 
sal (ow. * pv ) 


where ris the mass mixture ratio my em 
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2. Momentum 
By Assumption 3, the only force acting on the free- 
stream flow is that due to the pressure gradient. If M is 
the momentum flux at flow area A, the change in momentum 


flux across pressure increment dp is 
dM = —Adp (3) 
The momentum flux can be written as the sum of the 
momentum fluxes cf the gas and liquid. Thus, 
X= mh, V, + mV: ny 


If the flow were allowed to continue at constant 
pressure, Ve and vo would become equal to each other at the 
mass-weighted mean velocity V. Since for this process, 


aM = 9, the value of V is given by 


(sity + mu)V = m,V,+ mV; (5) 


V+ 1V; a 
aoc ae (3) 


“nus, the momentum flux can be written 
uM = (m, + m)V 
Since m_ + Ty is constant, the change in momentum flux is 


dM = (sn, + mud ee 


Substituting Eqs. (3) and (2) into Eq. (3), iV can be written 


Paes ay [a Dre 


= ler Pee pV; (3) 
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The slip ratio is defined as 
= V/V = (V,—-Vi/V Se 


This equation can be combined with Eq. (6) to give V_ and 


V, in terms of J: 


L 
aN SS 
v,=(1+ raz) ¥=0 (11) 
=ij—- ——)\7 = bv 
vi ( re)? 2 C12) 


The gas density can be expressed as 


pe = W,p/RT, 


where a is the effective molecular weight of the gas mixture 
and R is the universal gas constant. Eq. (13) is the defini- 
tion of the effective molecular weight Wes which is the 


quantity that gives the actual gas density when substituted 


Substituting Zqs. (11) - (13) into =q. (3), the dif- 


ferential momentum equation is 


oe - 2 RT. : r 
wad = dv = - (Gs +7) o (14) 


The guantities a and b are slowly varying because s is 
typically only 0.1 to 9.3 and slowly varying. The quanti- 


Les: fy T 5 wy and p, are also slowly varying. Integrating 
hed 
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Eq. (14) over a pressure increment Ap, for which a, b, r, 


To Wg? and o; are constant to within the desired accuracy, 


the change in 7? is 


2 RT,,, ptap/2 | texAp 
als = l+te aw, | p — ap/2 - ES) (15) 
All quantities other than pressure can be taken outside the 


integral and evaluated at their mean values (dencted by 


sub 
script m) correspondin to the mid-interval ressure 0 


Thus, 
- 2 
a 1+ tft. 
ry] ae 
RT,, [" ? dp. rm [ * 
= pouige ape (16) 
*|\ aw, if . » boul, * 
i- 
3 


ap 
oar 
Pea - at (ATs, t\ 
. re l¢+r aw,p * bp)? 
mee} 


(17) 


inal form of the momentum eauation. 
energy 


The enthalpy change of the mixture between state 1 
(the beginning of pressure interval Ap) and state 


A 
“ 


(the end 


31 


of the interval) can be evaluated in two steps: (1) phase 
change at 9), Ty? le and (2) change to Pos T52) Tho? at 
fixed compositicn. 


The enthalpy change for step 1 is 


enthalpy required 


aH amount of A to vaporize and 
aS vaporized x heat unit mass of 
a: From: Gi. tec T 
L gL 
| 
amount of B enthalpy required 


+ vaporized X to vaporize and 
heat unit mass of 


3 from ae to. T 


=] $1 
Kinetic energy 
amount of A required to 
+ and 3 Xx accelerate unit 
masc from V, to 
V Tl 
By 
or 
4H, = (ta, = ta, ) [Le, = Cay, (T,, ~ T:,)] ‘ 
* (Mo,. i me, .) [Lo, ba a (T,, os T.,)] oe 


+ (ty, — 1%, )(Vi.— Vi )/2 


where L and c are latent heat and specific heat, respectively. 


Y 


Introducing mere compact notation. 


4H, = Arty, (Le, + ce, 8:7) : 
(19) 


+ Arty, (Ls, + 6, 3:7) + md 


a 


The enthalpy change for step 2 is evaluated from the 
temperature, pressure, and velocity changes, with properties 


avaluated at mean T and p for the interval. 


AH, = ty, [614 (To, — Ta) + (Vin — Vi] 


+ tu, [er (T:, — Ti) + Saas 


(20) 
+ 1 3 _ Vv? | 
3s (Vi, 1) 
av} 
= ti, (4, AT, + =") 
7 A av? 5 
Au, (61.47: + 22 + (21) 


3y Assumption 3, no work is done by the free-stream flow 


and no heat is transferred to it. Hence, 
4H, + 4H, =0 C2 


Substituting Eqs. (19) and (21) into Eq. (22) and 


solving for AT, gives the energy equation for the mixture: 


-_1f4% Ap | avi 
AT, = Se [mgt tt (cnt + SE 4 SE ) 


am,3,V: Ame, 
+e + So (4, + a, 87) 


am, 
+ =, (t, + Coy, a7) 


4. Drag 
Although no force other than pressure acts on ‘the 


Sree-stream fiow as a whole, a drag force exists between 
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the phases. Hence, a second momentum equation must be 


witten using as the ecntrol volume the boundary between 


the phases. 


The two forces acting on each liquid drop are the 


buoyancy due to the pressure 
the relative gas velocity. 
the mass ti 


es the acceleration of the drop. Thus, 


single drop. 


gradient and the drag due to 
The sum of these is equal to 


for a 


cynamic pressure of drag frontal area 
relative gas flow x coefficient X ef drop 
i volume XY pressure g mass XY acceleration 
of drop : gradient drop drop of drop 
or 
1 «D? wD* dp 
% pylVi\Ve) CoAG ~ Ede 
G ee (24) 
" dV 
= (Fe) E) 
The absolute value sign in the first term makes che 
drag force oositive when : ee y, and negative when Vs < T, : 
So ot -_ 
’ Solving Hq. (24) for avy > 
( 
{ « 
i 
Sp, |s| sV2Codx _ _¢P 
av, = Selle - or 
pivi 
eas (25) 
Zu 


Differentiating Eq. (12), qv, can also be expressed in 


terms of s, r, and ¥V. Thus, 


= v sdr a ds 
dV, = bdV + (aa I+r (28) 
Solving for ds, 
bend? | sde (Lt sav, 
alana came Foe (27) 


Substituting dv, from Eq. (25), noting that dV = 


Pc ae : ee bos sp ate : 
dv/2V, using Eq. (12), and writing for a finite increment, 


results in, 
oe ba ce! + fm) av? + (l+ fu) Ap Saf 
2Vv2, bapi,V2i +t. 


$p,,, [su $aCo,, (1 = Ta) ax 
= 4 impr, D 


As 


This is the drag equation employed when x is specified 
as a function of oa. 
Solving Eq. (23) for A x yields the required alter- 


native equation: 


os 4D ee 53 pr, av? 
3p. | | SaCrV? [ SE Haan 
be pi, V?, SmAT eo 
l+ te fe ~*)| 
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S. Heat Transfer 

Although no heat is transferred to the mixture as a 
whole, heat transfer exists between the phases. Hence, a 
second energy equation must be written using as the control 
volume the boundary between the phases. 

The work dW done onthe liquid is that due to drag 
by the gas. (Only work done by shear or shaft forces is 
included in dW when writing the First Law for a control 
volume). Multiplying Eq. (24) by the number flow rate of 
drops N = bi, /mD°o, the drag force F, on that quantity of 
liquid is 


N m, d : dv 
Fe = Pa |Vs|V.CowD! = Ft SP + uV, (30) 


The work done on the liguid is 


(40 4 i 
~—dW = Fx =m ra 7) 


The heat dQ transferred from the liquid is made up of 
two parts: (1) the convective cooling due to the temperature 
difference between the liquid and gas and (2) the evaporative 
cooling due to the latent heat supplied to the liquid vaporized. 
The convective cooling is 


-dQ. = hAWN(T, ~ Tad 
(32) 


36 


where h is the heat-transfer coefficient, Aa 7 nD* is the 


surface area of a drop, and dt = dx/V, is the time required 


to traverse dx. Thus, 


6hrii(T, — T.)dz oe 


meee Dp.V, 


The evaporative cooling is 


—dQ, = Lidm,, + Lidm,, eR) 


The change is enthalpy of the liquid over the pressure 
increment dp is 


=a ar, + 22 4 SH) 
iin) ne Tear (35) 


into the 


a 


Substituting Eqs...€3l)s,. (33),.134)5 and: C35) 


steady-~flow energy equation dQ - dW = dH, the result is 


6hin8Tdz a ohn 
De - L.dm,, = Lidm, = micidT, 


where §T = T_ - Th 


Nrit: : for a finite interval, the final form of the 


fie 


drop heat-+transfer equation is 


-_ 1 6h8..T ax pay Amey or, am, 
aa i Cte [ DpiV ie . m1, i May ( 37) 
uh 
é al 


Equations (1), (17), (23), (23), and (37) are the five 


equations that must be solved simultaneously to obtain the 
values of the five dependent variables D, ee Ths Le and 

vr as a function of the independent variable p. To carry out 
the solution all quantities in the equations must be expressed 


in terms of these six variables. 


C. WALL SHEAR AND BOUNDARY LAYER 


For a two-phase nozzle, the momentum flux of the friction- 


less nozzle flow is that given by 


Yois 
m, Pr 
foo = 
p= AV OTt+h 
(38) 
where r_is the ratio of gas flow area to licuid flow area 


of 3 at a station where the nozzle wall radius is a is 


where M,. is the momentum flux of the real Flow with friction. 
4 


The skin-friction coefficient can be defined using the same 


quantities as single-phase flow 


C= Qe 
1 Oe (40) 


where T, is the wall shear. It will be shown that a valid 


Ce value can be provided. 


The buundary-layer momentum equation then becomes 


(41) 
1+ (8/8) = 1 = 1 
~— 6 | ——ee—— dV + —=d(pV) + —dR,, 
[ . + pale) +e | 
| where 6* is the displacement thickness, i.e., the distance 
the wall must be moved outward to give the same flow rate 
as with frictionless flow. 
1 
Assuming a 5 power velocity profile and no density vari- 


ation, the shape factor $*/3s is obtained from 


: 5* [ * — (y/sy"} dy 
7 -— 


| ° (y/8)" [1 — (y/8)") dy 


=i} 


where 3 is the velocity thickness of the boundary layer. 


+> 


Noting that dV can be written av*s2v, and that 4(9'T) = 


a(m,/A), the finite-difference form of Za. (41) is 


Wall shear in homogeneous two-phase flow has been found 


to be egual to thet which would exist with pure Liquid at 


equal velocity and boundary-layer thickness, multiplied by 


the wetted wall fraction: 


Cy, 2 Ay C. p V: 
FOV = TET (42) 


is the skin friction coefficient for liquid at 4 


where a 
Reynolds number of 
Rs = pivis 
ry 


For a 1/7-power profile, the velocity thickness, 


A convenient relation for Ce as a function of Re is the 


Shultz-Grunow relation which can be written 


a 8 ee 
Cy, aE (logio Rs ~” 0.425)" 5° 


pee 


Comparison of Eqs. (40) and (42) shows that C, can be 


written 


Thus, the final form of the boundary-layer momentum 


equation 8 


fmm Guar ~ 6, (Sp. - SA + Be) 
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Let Vs be the mean velocity of the flow including the 


boundary layer. Then, from Eq. (39), 


= 


m,Vs = M,= mV — Ixy. p'V26 


Rence, employing Eq. (38), the mean exit velocity includ- 


ing the boundary layer is 


By the definition of the displacement thickness, the flow 


rate is reduced by the throat displacement thickness 5¢ to 


( 287 
» =m i jl-— 
™s mr Yo, 


D. NOZZLE THEORY SUMMARY 

The preceding equations form the basis for the mathematical 
model which is used to predict, based on inlet conditions, the 
exit velocity, and temperature of the mixture. These equations 
also form the basis for the model which vrovides the optimum 
nozzle shape given a set of inlet conditions. Some additional 
relationships are, however, reguired. These are: 

1. Phase properties - te establish the mass ratio, mass 

Flow rate ratio of gas to liquid, and the thermal 


conductivity of the mixture. 


4 
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2. Liquid drop drag coefficients. 
3. The liquid drop heat transfer coefficients. 
4. S3oundary layer momentum thickness and displacement 
thickness. 
5. Skin friction coefficient. 
These five additional relationships are developed in detail 


in Reference [1]. 
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IIIT. COMPUTER PROGRAM DUAL-PHASE NOZZLE 


The computer program employed in this study is based on 
a program developed by Dr. G. Elliott of the Jet Propulsion 
Laboratory in Pasadena, California. The program was updated 
and converted for use on the Naval Postgraduate School 
computer. The Dual-Phase Two-Component program employs the 
theory in Section II. The program is written in Fortran 
computer language and can be compiled using a Watfiv or For- 
tran IV compiler, 

This program has been utilized in dual-shase nozzle 
analysis and to previde values for comparison with ‘he ex- 
perimental results. To use the computer program the inlet 
conditions have to be specified. The flow conditions are: 
inlet pressure; mass ratio; inlet temperature of the gas and 
liquid; inlet velocity of the gas and liquid; ‘total mass flow 


a 


rate; and nozzle exit pressure. Section III 8, shows speci- 
fied detaiis for data input. 

There are two options that can be chosen. The first is 
prescribed pressure-versus-distance option MOP=0. The ores- 
sure profile ?(X) is selected corresponding to the adopted 
nozzle contour. If the pressure-versus-distance is used a 


P(X) input table is required. This profile is develcped 


from the actual measured pressure values in the experiments. 
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(See Appendix A for sample program.) The second option con- 
sist of an optimum nozzle contour option MOP=1. This option 
is useful only when the liquid drop diameter is constant. 

The dual-phase two-component computer program is a struc- 
tured program with thirteen subroutines controlled by a main 
program. This arrangement improved the programming process 
through better organization and programming notation. 

The control point of the dual-phase two-component com- 
puter program is the "main section." It controls the flow 
path and operation of all input data, property tables, and 
calculations. It accomplishes this by calling the thirteen 
subroutines at the appropriate times, saving wanted data 
in files, and printing out desired information. 

One of the most important subroutines which inputs in- 
formation is the "INTRP” subroutines. INTRP controls the 


property table inputs. It reads in four two-dimensional 


ct 


ables and fourteen one~dimensional tables. These inputs 


t are the properties of the gas and liquid phases of Doth 


components of the flow in the nozzle. The subroutine writes 
the values of these tables into a file and retrieves appro- 
priate values from that file. INTRP can also interpolate 
for values used throughcut the entire program. 

Input of case data is controlled by subroutine "Sect 1.” 
Identification information and case heading information is 
read and printed for each instance. Sect 1 also places the 


pressure vs. distance profile, if specified, in an array. 


4u 


"Sect 2" through "Sect 6" are the subroutines which cal- 
culates the flow data. Sect 2 sets the initial conditions 
indicated in the input, and begins the iterations. Sect 3 
computes initial flow rates of both components; the initial 
area of the nozzle; slip velocity; mean free stream velocity; 
and slip friction. Sect 4 computes the changes in flow par- 
ameters and new distances. It then begins to calculate new 
conditions such as flow rate, temperature, velocities, sur- 
face tension and mean area. Sect 5 is the binary cut con- 
vergence routine and computes mean boundary-layer parameters. 

If a problem is diagnosed in any subroutines and "diagno" 
is called, it will p: -nt all output parameters calculated to 
that point. It also does the same if there is a convergence 
problem. 

There are two subroutines that output calculated data, 
subroutine "Write"™ and subrcutine "Outzut." Subroutine 
"Nrite" will send output information tc the printer for a 
hard paper copy. Subreutine "Output" reads and stores the 
Gutoat-6n.a file, 

The two-chase two-component program has been written 
with comment statements in the text of the crogram. These 


will allow for a more understandable and, therefore, a more 


wn 
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easily modified program. For specific detail 
of these subroutines, see Appendix M. 
Tne dual-phase iozzle program was tested for correct 


output. Sample data and results were obtained from 
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a~--—ea b heme oe on wo mn ne ee eee 


— 


Dr. David Elliott, were inputted into the program, and 
executed. The output was compared to the sample data. The 
program produced duplicate results. 

The program begins by storing fluid property tables and 
reading in all input data. All nozzle inlet conditions are 
computed. The program then proceedshalf a pressure step at 
a time. At the middle of each pressure interval, the changes 
in quantities across the interval are computed using the pro- 
perties interpolated from the table for that pressure, and 


for the existing temperature. The change in slip is found 


th 


if the pressure profile P(X) is specified. At the end of 
each pressure step, the flow conditions are updated and 
initial conditions are determined for the next step. The 
dropsize is reduced at the point when the Weber number ex- 
ceeds six. The flow conditions are printed if the pressure 
is one selected for output. The computation continues until 
the last pressure step has been completed and flow conditions 


at the smallest flow area encountered are printed as the 


threat conditions. 


A. PROPERTY TABLES 


2. Heat capacity of component "A" vapor in 3TU/L3M-=R 
is a function of temperature and pressure. The two-dimen- 
Siomal tables are entered row-wise. At least two cards are 


necessary to specify a row and at least two rows must be 


entered, 


me ‘sos inn nen nere 


Card 1: (format 6E£12.6) 
cols. 
1-12 temperature (R) 
13-2421.0 if this is the last temp for 
this table 
Card 2: (format 6£12.6) 
cols. 
1-12 pressure (psi) 
13-24 heat capacity (BTU/LBM-R) 
25-36 pressure (psi) 
37-48 heat capacity (BTU/LBM-R) 
49-60 pressure (psi) 
61-72 heat capacity (BTU/LBM-R) 

The maximum entries of temperature are 35 values. For 
each value of temperature, the maximum number of entries of 
pressure and heat capacities are 35 values. Each row of this 
table will be terminated with the pressure and heat capacity 
equal to to”. These two values are not counted in the maxi- 
mum of 35 entries/rows allowed. 

The program shown in Appendix 8 can be used to determine 
the values of heat capacities. The program structures its 
output in the format needed for the table input. It uses 
input data obtained from Reference 3. The input data must 
be placed in a two-dimensional table. This table is used 
in the program to interpolate the values needed for output. 


input data must be formatted as follows: 
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Card 1-16 (format 10F7.4) 
cols. 
1-7 temperature (R) 
at this temp the following is entered: 
8-14 heat capacity at .01 P 
15-21 heat capacity at .4 Pa 
22-28 heat capacity at .7 Pa 
29-35 heat capacity at 1.0 Pa 
34-42 heat capacity at 4.0 Pa 
43-49 heat capacity at 7.0 Pa 
50-56 heat capacity at 10.0 Pa 
97-63 heat capacity at 40.0 Pa 
64-70 heat capacity at 70.0 Pa 
The temperature must be entered with 
increasing value. 

2. Heat capacity of component "B" gas, BTU/LBM-R is a 
function of temperature and pressure. This two-dimensional 
table has the same format as part Al above. 

3. Molecular weight of component "A" vapor is a funce 
tion of temperature and pressure. The two-dimensional 
tables are entered row-wise. At least two cards are neces- 
Sary to specify a row and at least two rows must be entered, 

Card 1: (format 6£12.6) 
cols. 


1-12 temperature (R) 
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13-24=1.0 if this is the last temp for 
this table 
Card 2: (format 6E12.6) 
cols. 
1-12 pressure (psi) 
13-24 molecular weight 
25-36 pressure (psi) 
37-48 molecular weight 
49-60 pressure (psi) 
61-72 molecular weight 

The maximum number entries of temperature are 35 values. 
For each value of temperature, the maximum entries of pres- 
sure and molecular weight are 35 values. Each row of this 
table will be terminated with the pressure and molecular 
weight equal to 10°, These two values are not counted in 
the maximum of 35 entries/rows allowed. 

The program shown in Appendix C can be used to deter- 
Mine the values of molecular weight. The program formats 
its output in the format needed for the table input of 
the two-component two-phase computer program. It uses 
input data obtained from Reference 3. The input data 
must be placed in a two-dimensional table. This table 
is used in the program to interpolate the values needed 
for output. The program used in Figure 8 can only be used 


with ideal gases. Input data must be formatted as follows: 
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Card 1-16 (format 10F7.4) 
cols. 
1-7 temperature (R) 
at this temp the following is entered: 
8-14 density at .01 P 
15-21 density at .4 Pa 
22-28 density at .7 Pa 
29~35 density at 1.0 Pa 
34-42 density at 4.0 Pa 
43-49 density at 7.0 Pa 
S0~56 density at 10.1 Pa 
$7-63 density at 40.0 Pa 
64-70 density at 70.0 Pa 
The temperatures must be entered with 
increasing value. 

4. Molecular weight of component "B" gas is a function 
of temperature and pressure. This two-dimensional table 
has the same format as part A3 above. 

5. There are fourteen one-dimensional tables. The 
one-dimensional tables are entered in the following format 

for Z(T)): 
Card 1: (format 6£12.6) 
cols. 


1-12 Ti-2, oR 


25-36 Ti-1l 2 <i «< 50 

37248 Zi-01 

49-60 Ti 

61-72 Zi and Tj-1 < Tj for 2 <j < 50 


Each table is terminated by two consecutive entries of 
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10°, i.e., Tk = Zk = 10° (1.0E5 right adjusted in the field). 


The fourteen one-dimensional tables are (in order requested): 


nee CAL(T) heat capacity of component A liquid, 
BTU/L3M OR 
2. CBL(T) heat capacity of component B liquid, 


BTU/LBM oR 

cn LA(T) latent heat of vaporization for com- 
ponent A, BTU/LBM 

4, L&(T) latent heat of vaporization for com- 


ponent B, BTU/LBM 


Drs PBOC(T) vapor pressure of component 8 

6. ROAL(T) density of liquid component A, 
LBM/FT3 

7. ROBL(T) density of liquid component 8B, 
LBM/FT3 

3. KAG(T) thermal conductivity of component 


A gas, BTU/FT HR oR 
9: KBG(T) thermal conductivity of component 


B gas, 3TU/FT HR oR 


Sete, F = em eee ee re oe 


10. 


12. 


14. 


VIAL(T) viscosity of liquid component A, 
LBM/FT HR 

VIBL(T) viscosity of liquid component 8, 
LBM/FT HR 

VIAG(T) viscosity of gas component A, 
LBM/FT HR 

VIBG(T) viscosity of gas component 8, 
LBM/FT HR 

SIG(T) surface tension of liquid component 


B, DYNE/CM 


sn 


Appendix D gives sample property table for input. 


3. CASE INPUT 


A blank card must separate the property table from the 


data set decks following. 


Card 1: (format 4A4, A2, 3A4, 7A4, A32, 215) 


cols. 
1-132 date 
19-30 case number, mav te any 
alphanumerie fata 


32-60 identification 


Oo 


1-66 NS, number of vressure steps per 
printout 

(right justified integer).* Use -1 if new 
property tabies follows. 

537-72 NP, number of printouts (right 


justified) 
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Card 2: (All integers right justified in field) 
format 1116) 
cols. 
1-8 MBU, =0 constant droplet size, =1l drop 
breakup 
7-12 MOP, =0 X(P) table to be supplied, =1 
X(P) is to be optimized 
13-18 MGEOQ, =0 circular, =1 annular 
19-24 NDS, maximum number of S iterations 
25-30 NSO, maximum number of So iterations 
31-36 NB, maximum number of TB iterations 
37-42 NNS, first setting of step counter 
Card 3: (format 55F12.5) 
cols. 
1-12 DP, pressure step size, negative for 
decreasing, psi 
13-24 RC, mass flow ratio 
25-36 PHI, angle of annular nozzle axis, deg 
37-438 RAXO, annular nozzle axis, in. 
49-60 EMT, total flow rate *P final = PO + 
(NS #NP-NNS-1)*DP4+DP1 
Card 4: (format 3£12.5) 
cols. 


1-12 H, inverse Henry's Law ccenstant, psi 


3-24 ALAM, Lagrangian multiplier 
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25-26 DPL, first pressure step size 
37-48 WAL, molecular weight of liquid a 
49-50 WBL, molecular weight of liquid b 
61-72 SA, Sutherland constant for component 
A, OR 
Card 5: (format 6£12.6) 
cols. 
1-12 SB, Sutherland constant for component 
B, oR 
13-24 DO, initial drop diameter, in. 
25-26 PO, initial pressure, psia 
37-438 TGO, initial gas temperature, OR 
49~60 TLO, initial liquid tempeature, oR 
61-72 VGO, initial gas velocity, FT/S 
Card 6: (format 6512.6) 
cols. 
1-12 VLO, Initial liquid velocity, FT/S 
13-24 THOO, initial momentum thickness cf 
outer wall boundary layer, in. 
25-36 THIO, initial momentum thickness of 
inner wail boundary layer, in. 
37-48 EDS, convergence criterion for $ 
49-60 ESO, convergence criterion for So 
§1-72 EB, convergence criterion for T, oR 
If MOP=0, the following cards are present: 


in 10 Card 7: (format 7A4, A2) cols. 
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1-30 X(P) table identification (and 
alphanumeric data). 
Card 8: (format 6£12.6) 
eols. 
1-12 pressure, pi-2, psia 
13-24 distance, xi-2, in. 
25-36 pressure, pi-1l G8 2.6.79 
37-48 distance, xi-1 
49-60 pressure, pi 
61-72 distance, xi 
The last two entries are 1.0£5 and 1.0E5 right adjusted 
in their fields. The table must be monotonic increasing or 
decreasing. New property tables may be used by putting -1l 
in cols. 61-66 of Card 1, and following this with new 
property tables and data sets. Appendix E is a sample 


input data. 


C. OUTPUT 
For each case, the case identification is printed followed 
by the input parameters. If MOP=0, the X(P) table forms a 


part of this output. The following output then appears. 


Lp. ox distance, in. 

Ze -P pressure, psia 

3. R mass flow ratio 

4. vb mean free-stream velocity, ft/s 
5S. a flow area, in. 2 
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6. tb gas temperature, or 

7. ti liquid temperature, or 

8. vg gas velocity, ft/s 

9. vi Liquid velocity, ft/s 

10. vs slip velocity vg-vl, ft/s 
| is eae slip fraction vs/vb 
| 12. d drop diameter, in. | 
| 13. rv ratio of gas volume flow to liquid volume flow 

14. ra ratio of gas flow area to liquid flow area 


15. alpha mass fraction of component a dissolved in 
liquid 


16. beta mass fraction of component b vapor in gas 


17. mg gas flow rate, lbm/s 
13. mil liquid flow rate, lbm/s 
19. rog gas density, lbm/ft3 
20. rol liquid density, lbm/ft3 
21. wag molecular weight of component a gas 
t 22. whg molecular weight of component b gas 
| 23. wg mean molecular weight of gas 
24. pa partial pressure of component a, pDsia 
25. pb partial pressure of component b, psia 
25. la latent heat of vaporization of component a, 
btu/1lbm 
27. 1d latent heat of vaporization of component b, 
tu/lbm 
56 
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28. 


29. 


sigma liquid surface tension, dyne/cm 


egm specific heat of gas (at midpoint of pressure 


step), btu/lbm of 


30. clm specific heat of liquid, btu/lbm of 


31. vigm viscosity of liquid, lbm/ft hr 


32. vilm viscosity of gas, lbm/ft hr 


33. kgm thermal conductivity of gas, btu/hr ft of 
34. rem reynolds number of flow over drops 
35. edm drag coefficient of drops 
36. hm heat transfer coefficient of drops, 
btu/hr ft2 
37. rb mass flow ratio after velocity and temperature 
equalization 
38. ab flow area after equalization, in. 2 
39. tb temperature after equalization, or 
4Q. rvb volume flow ratio after equalization 
41. alphb alpha after equalization 
42. betab beta after equalization 
43. mgb gas flow rate after equalization, lbm/s 
44, mlb liquid flow rate after equalization, l1bm/s 
45. vilb liquid viscosity after equalization, 
lbm/ft hr 
46. rogb gas density after equalization, lbm/ft3 
47. rolb liquid density after equalization, lbm/fts 
48. wagb molecular weight of component a gas after 
equalization 
57 
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43. wbgb molecular weight of component b gas after 


equalization 
50. wgb mean molecular weight of gas after equalization 
S1. pab partial pressure of component a gas after 


equalization, psia 
52. pbb partial pressure of component b gas after 


equalization, psia 


Ceres separator friction parameter 

54. ref separator film reynolds number 

55. yo distance from nozzle axis to outer wall, in. 
56. wom angle of outer wall relative to axis, deg 

57. tho momentum thickness of outer boundary layer, in. 


53. delo velocity thickness of outer boundary layer, in. 
59. delso displacement thickness of outer boundary layer, 
in. 
60. redom reynolds number of outer boundary layer 
51. ecfom skin friction coefficient of outer boundary 
! layer 


52. twom shear stress on outer wall, psi 


63. vbd mean velocity including boundary layer, ft/s 
34. wim angie of inner wall relative to axis, deg 
30%, tHL momentum thickness of inner boundary layer, in. 


66. deli velocity thickness of inner boundary layer, in. 
57. delsi displacement thickness of inner boundary 


layer, in. 


68. redim reynolds number of inner boundary layer 
69, cfim skin friction coefficient of inner boundary 


| layer 


70. twin shear stress on inner wall, psi 


Jil. nna number of iterations required to optimize X(P) 
72. nis number of iterations required to converge cn 
s or so 
Tae RO number of iterations required to converge on 
tb 
$3 
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Iv. PERIMENTAL SYSTEM 


The experimental system can be grouped into three sub- 
systems. These are: 

a) nozzle 

b) air supply system 

¢) liquid injection system 
Each subsystem is described in the following sections. 


Figure 11 is an overall system schematic. 


A. NOZZLE 

The nozzle has a convergent-divergent flow passage. It 
is 12 inches long with a variable exit area. The exit area 
can be varizd from .45313 square inches to .84375 square 
inches. The pivot point is located 1 inch above the throat. 


This causes the throat to vary when the exit is varied. 


Since the change in the throat is negligible, it will be con- 


idered to be constant. It has a throat area of .45 square 


inches. The inlet area is 1.625 square inches. The throat 


is located 4 inches from the inlet. The nozzle is constructed 


by sandwiching two 1/2 inche thick machined aluminum nozzle 
profile plates between 1/2 inch plexiglas plates (Fig. 12). 
The aluminum nozzle plates are located at the end of a 30 
inch long entry section and are easily adjustable. Figure 


13 shows a close up of the aluminum section of the nozzle. 
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The inlet section provides the necessary space for pressure 
taps, liquid injection tubes, and other instrumentation and 


configuration options. 


B. AIR SYSTEM 

The air system is shown in Figure 14. Two 117 cubic 
feet tanks supplied with compressed air from an Ingersol 
Rand two-stage air compressor provides the air storage vol- 
ume required to support the nozzle's operation. Each stor- 
age tank has a pneumatically control Norgren gate valve at 
its exit. These valves are opevated by a nitrogen actuator 
and controlled by a pressure regulator. The nitrogen is 
regulated to 40 psi control pressure which will cpen the 
Norgren valve. 

The nozzle is supplied with air via a 3" i.d. pipe. The Ls 
air supply to the nozzle is controlled by a solenoid actuated 
nitrogen operated 3 inch ball valve. The nitrogen is sup- - 
plied via a regulator. By varying the nitrogen supply pres- 
sure to the ball valve, supply air pressure to the nozzle 


can be controlled. Air flow to the nozzle is measured with 


a standard ASME orifice plate. Figure 15 shows the dimen- 
sions of the orifice. The orifice is a model D-19512 with 


a 9.929 inch bore. tO 


C. LIQUID INJECTION SYSTEM €] 
Tne liquid iniection system, Figure 16, is supplied by 


house water and is further pressurized Sy an Aurora electric 
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pump model/type 317696. The pump is rated at 50 gallons per 
min. The pump supplies pressurized water to the injection 
tube via two flowmeters (rotometers). The flowmeters are 
F&P co. precision bore flowrator tubes. One rotometer is 
rated at 1 to 12 gallons per min. and the other at .6 gallons 
per min. 

The liquid injector is a 0.25 inch brass tube inserted 
in the 3" i.d. air supply pipe just upstream of the flange 
connection to the test section. The intector tube is drilled 
with sixteen 1/16th inch diameter holes facing the test sec- 
tion entrance. The drilled holes were made as small as 
possibile consistent with achieving a significant Liquid 


mass flow rate. 
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V. INSTRUMENTATION SYSTEM 


The instrumentation system is designed for automated data 
acquisition, analysis, and display. A schematic is depicted 
in Figure 17. Where ever possible, test operations and 
sequencing are under direct computer control. Each parameter 
measured involves an appropriate transducer, excitation source, 
and calibration procedure. The maior instrumentation sub- 
systems are: 

(A) Pressure Measuring Transducers 

(B) Nozzle Thrust Force Block 

(C) Flow Measurement Devices 


(D) Data Acquisition/Control System 


A. PRESSURE MEASURING TRANSDUCERS 

Pressure measurements are made in fourteen locations 
throughout the experimental apparatus. Eleven Micro Twitch 
14O0PC pressure transducers model PK 87633 are placed on the 
Nozzle assembly to measure pressure at various axial posi- 
tions. ‘iIpecifications for this model transducer are shown in 
Figure 18. The first pressure tap is located at one half 
inch from the inlet along the axis of the nozzle. The 
remainder are placed at one inch intervals toward the 


nozzle exit. These pressure taps are connected to the 


pressure transducers via a 1/4" o.¢. plastic tubing. The 
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pressure transducers are mounted on the side panel of the 
nozzle. See Figure 19 for locations. Each transducer has 
an identification letter corresponding to its position. 
Table I is the ID code for the transducers. 

At the nozzle inlet a Micro Switch 200PC pressure trans- 
ducer model PK 87690 is placedin the inlet line to measure 
the dual-phase mixture inlet pressure. Specifications for 
this model transducers are shown in Figure 19. Gas supply 
pressure is measured using a 200PC pressure transducer 
located at the orifice outlet. This transducer is also 
used to help measure differential pressure across the ori- 
fice. Another 200PC transducer is placed at the inlet of 
the orifice. The data acquisition program calculates the 
orifice differential pressure and then the air mass rate 
using the output from the above two transducers. 

The pressure transducers provide output voltage propor- 


tional to applied pressure. These operate from a single, 


positive supply voltage ranging from 8 to 20 VDC. The 
supply voltage was maintained at 15VDC. 

Each pressure transducer sends back a signal proportional 
to the input pressure. The signal is converted in the A/D 
converter to a digital signal which can be read by the HP9826. 
The pressure transducers were calibrated using a known pres- 
sure source. Appendix F depicts a sample calibration pro- 


gram written for a HP9826 computer. This program reads ten 
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208PC SPECIFICATIONS of 6.0 2 0.81VOC, 23°C (uniese otherwise noted) 
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Table I. ID Code for Instrumentation 


Letter Designation Description 
A 10.5" from inlet 
B 9.5" from inlet 
Cc 8.5" from inlet 
D 7.5" from inlet 
E 6.5" from inlet 
F 5.5" from inlet 
G 4,5" from inlet 
H 3.5" from inlet 
I 2.5" from inlet 
J 1.5" from inlet 
K 0.5" from inlet 
L inlet nozzle pressure 
M force-block signal 
N force-block excitation voltage 
0 orifice exit pressure 
Pp orifice inlet pressure 
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Pressure = A + B(volts) + ¢(volts)* + O(volts)? 


pressure 
transducer A 8 C i} 


ean 


A -6.069 3.926  -.17336 ** .01123 
B 7.059 4.234  -.22921 = 0415 
c -.6491 3.934 -.16025 .00945 
D “7.056 4.692 -.34118 01896 | 
E 27.009 4.152 = -.20717— «.01237 | 
F -6.389 3.933  -.16055 .00957 

6 “6.589 3.956  -.16427 .00967 

Xt 6.588 3.908 -.15806 00942 a 
1 “6.595 3.949  -.16270 ~~ .00957 . 

F) -6.353 3.927 -.J6509  .00998 

K -6.203 . 3.846  -.14697  .00880 

L 20.207. 11.157 =.08591 —-_.00492 

0 20.668 11.068 -.07562 .01118 

P 18.827. 11.928 = =.17715— 01101 


Figure 20. List of Polynomial Coefficients 
for Pressure Transducer 
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values of pressure for each 140PC transducer and gives the 
mean and standard deviation for each from 0 to 45 psi in 5 
psi increments. Appendix G depicts sample output of the 
program for a given pressure. Appendix H illustrates the 
program used to calibrate the 200PC pressure transducers. 
The program works as above, except values are taken from 0 
to 60 psi at S psi increments. 

Data obtained during calibration source pressure is 
plotted for each pressure transducer (Appendix I). Each 
plot was curve fitted with a third order polynomial. Figure 
20 shows coefficients of the polynomial. These polynomials 
are used in the data acquisition/control program to convert 


transducer readings to pressure readings. 


B. NOZZLE THRUST FORCE-BLOCK 
The thrust produced by the nozzle was employed to deduce 
the nozzle exit velocities. The thrust was measured by in- 
strumenting a target plate in the exit flow field. Appro- 
priate screens were installed to prevent liquid “bounce 
' back." The jet momentum force on the target is acquired by 
a balance beam system shown in Figure 21. A Kistler-Morse 
force block provides an analog signal proportional to the 
nozzle jet momentum. 
The calibration of the force-block was completed by 


placing known weights on the force-block side of the balance 
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beam, (Figure 21), and recording the voltage produced. 
Appendix J shows results of this calibration. These results 
were curve fitted with a third order polynomial. The follow- 
ing is the volt-to-force conversion polynomial: 

FB = (1/453.6)*(-234 396448715 .28*VOLTS-11261 .54*VOLTS+33899 . 3*VOLTS) 

FN = FB*®LB/LN 

FN = FORCE AT NOZZLE EXIT (LBF) 

FB = FORCE AT FORCE BLOCK (LBS) 

LB = LENGTH FROM PIVOT TO FORCE BLOCKCIN.) 


LN = LENGTH FROM PIVOT TO DIRECTIONAL OBSTRUCTIONCIN.) 


C. FLOW MEASUREMENT 

1. Air mass flow measurement for the system achieved 
by measuring the inlet and outlet pressure on a D-10512 ori- 
fice with a 0.92 inch bore. for details on orifice see 
Section III. Air flow is calculated in the data acqui- 
sition program. The program uses the following equation 


obtained from References 4% and §. 


. = : Pi- I 2 
Maip = K Ay YV%S, po, CPLoPe) her 
K = CE 
hae tle 
B= 4/D 
Cs 0.80 
dy = AREA OF ORIFICE = nD? FT? 
ta 
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e. = 32.2 FT/SEC2 

P, = DENSITY LBM/FT3 

Pl = PRESSURE AT INLET IN LBF/FT2 
P2 = PRESSURE AT LUTLET IN LBF/FT2 
Moir? AIR MASS FLOW RATE LBM/SEC 


The discharge coefficient C is the factor that accounts 
for losses through the orifice. Since the values of C varies 
from .52 to .69 ror R, number from ao to 107 with B = .32, 

C will be considered constant. 

2. Water flow measurement was made using two rotometers. 
Calibration of the rotometers were made by measuring the 
time for a given weight of fluid flow. The mass flow rate 
was calculated and plotted versus the rotometer reading. 
Appendix K is the plot of the results. The plot was curve 
Fitted with a third order polynomial. The following is 
that polynomial: 

MH20 = ~ .0063268+.902097278*#RR~.00000658*RR-+1,1X10#RR° 


where RR = rotometer reading 


D. DATA ACQUISITION/CONTROL SYSTEM 

The heart of the data acquisition system is a HP9826 
small computer. The HP9825 communicates via a Hewlett 
Packard 3497A data acquisition/control system. This system 
gathers data from the pressure transducers and nozzle thrust 
force-block. It converts the analog signal to digital data, 
and stores the data in memory. Figure 22 shows a pressure 


cransducer to A/D converter channel connection, and Figure 23 


5 PTB ae rtqn et ise. se ewe Ss Pee 


uoT OoUUOD TeuUPYQ-o}-daonpsuPaL *2Z Panky 


T3ANNVH) N3LY3ANOD C/V ' 


603 [0A 04303 ,9x9 
aatzebeu 


9622 10A UOZ2eZ4,9xa 
3A$34S0d 


_ - punoib 


UINGSNVUL FUNSSIUd 


UOT}POSaUUOD Toauueyg-01-xo0TgG-a0u0g 


TANNVHD YW3LYIANOD O/¥ 


NOILVLIOX3 NOILVLIIXA 
JAILIOIN SAILISOd 


waa... 


49078 39u04 


"eZ aunsty 


MOT LNT 


ewe 


Seceld 


ca mn ee ee 


ee ee 


SUOTJOSUUOD Tayteauoy) G/y-0}-UOTJOaUULD sng 


"he suNnsTy 


UZLYIANOD O/V OL NOILIINNOD SNE ' 


21 ‘T3NNVH Ut ‘TaNNVHD OL T3NNVHD € ‘TaNNVHD 2 T3NNVHD =. = L T3NNYHD 


N° d 0 . F | 3 a 


Sl T3NNVHD YL TINNY €L T3NNVHD 9 TANNVHD § T3NNVHD ¥ ‘TINNVHD 


JOVLIOA 
NOILVLIOX3 I H 9 


e 
SU TANNYHD CL TANNVHD 9. ‘TaNNVH 6 TANNVHD 8 ‘T3aNNVH < ‘TaNNVHD 


a a SE AS, tS 


rat ‘ Woe os, BHye die. 


‘ Arad 
r a aed 


ne ne em ee ee 


a 


the force-block to A/D converter channel connection. Figure 
24 depicts the current connection to the data acquisition/ 
control system bus from the sensors. 

The data acquisition and analysis program is a fully 
automated system which integrates the use of the HP9826, 
HP3497A, and all sensor devices. The program can be divided 
into three sections: 

1. Data Acquisition - Data is acquired from 15 differ- 
ent sensor devices. All system parameters are ini- 
tialized and configuration types are inputted into 
the computer memory for future use. Ten readings 
are made om each sensor and a mean value for each 
is calculated and stored in memory. 

2. Data Conversion ~ Input data conversion is completed 
by using calibration information obtained on each 
pressure transducer and force~block. The program 
takes the stored mean values from memory, and 
converts the input voltage to a pressure or force 
using the third order calibration polynomial for 


each sensor. 


3. Data Analysis - This is the heart of the program 
in that all performance parameters are calculated 
and printed. It uses all of the above information 
to calculate mass flow rate of water, mass flow 


rate of air, total mass flow rate, thrust, mixture 


ratio, exit velocity and all pressure informa- 

tion. The following equations are also used: 
Exit Velocity = (Thrust/Total Mass Flow Rate)* g, (ft/sec) 
Mixture Ratio = Mater’ Mair 


+ m 


Total Mass Flow Rate = MN vetan Soe Sooo 


Appendix L is the data acquisition program and 


sample output. 
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VI. EXPERIMENTAL RESULTS 


Six test runs were made using the equipment described in 
the previous sections. It was found that the inlet pressure 
to the nozzle remained nearly constant for each setting. 

The testing was conducted by maintaining a constant air 
control valve setting and varying the mass flow rate of 
water, The regulator outlet pressure of the nitrogen bottle 
which controls the air valve was set at 30 psi and 15 psi. 
These settings correspond to nozzle inlet pressure of 35 + 
1.0 psi and 29 + 1.0 psi respectively. For each nozzle in- 
let pressure, experimental data was obtained at varying 
mass ratios of water-to-air in the range of 2 to 13. These 
experiments were conducted with three discharge area ratios 
(Ag niet/Aout?: Those ratios are: 3.586, 2.600, and 
1.9259. The results of the six experimental tests are dis- 
played in Tables II through VII. 

Inlet air was limited to a maximum value of 40 psi. 

The 140PC pressure transducer's maximum output pressure is 
40 psi, (see information on 140PC transducer Figure 21). 
The experimental system has the ability to reach higher 
pressures. The present configuration may be overated to 


70 psig. 
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VII. DISCUSSION 


The results of this investigation are in two parts; the 
experimental test results and the computer outputs based on 
initial conditions similar to those of the experiments. 

The key variables are: the nozzle overall area ratio 
Ap, liquid/air mass flow ratio, nozzle exit velocity, nozzle 
Supply pressure, nozzle thrust, and the nozzle axial pres- 
Sure profiles. Figures 25 through 34 present the exit vel- 
ocity vs. mass ratio results for the experiment and computer 
output. Figure 29-34 illustrates the comparison of the ex- 
perimental results and the corresponding computer output. 


Figure 35 illustrates the variation of nozzle thrust with 


mass ratio for different area ratios and inlet pressures. e 
Figures 36 and 37 present the axial pressure profiles. 
In all cases the exit plane velocity decreases as the 
| liquid/air muss ratio is increased. In the low mass ratio 
range (i.e., less than ~ 5) the velocity decrease is very 
pronounced. Past a mass ratio value of ~ 10 the velocity 
of the mixture becomes relatively insensitive to the mixture 
ratio. This is as wouli be expected. For the same inlet 
conditions there is a f:.xed amount of energy available for 
conversion to its kinetic form. In the nozzle process the 


energy is conserved and thus an increasing mass ratio is 
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manifested in a decreasing exit plane velocity. In all cases 


the agreement between the experimental tests and the analytic 
predictions are with in 2 10%. 

There also was apparent and a consistent trend with re- 
spect to the nozzle area ratio and the exit velocity. As 
the nozzle exit area was increased the exit velocity of the 
mixture decreased. The effect seems to be more pronounced 
at the lower mass ratios. This trend was confirmed by ob- 
servations of the actual flow field in the nozzle. if the 
exit area was increased significantly (beyond the max area 
used in those tests) there was evident an abrupt and severe 
separation within the diverging portion of the nozzle 
passage. This was accompanied by a drastic decrease in 
exit velocity as evidenced by the output from the thrust 
target. It appears that the diverging portion of the 
nozzle, at a certain ~.int, starts to behave as a subsonic 
diffuser and hence experiences an adverse pressure gradient. 
This reasoning is in part confirmed by the pressure profiles 
identified in Figures 36 and 37. 

The relationship between measured thrust and mass ratio 
(Figure 35) indicates a slightly increasing trend. This may 
be explained by considering the following. Each set of data 
points (at constant inlet air pressure and constant nozzle 
exit area) is developed by varying the mass ratio. This, 
in turn, is achieved by increasing the liquid rate by in- 


creasing the liquid supply pressure. The net result is an 
Diy ¢ 


( 
ww 


increase in the liquid inlet velocity or an increase in the 
inlet energy level. Thus a particular data set is not truly 
at a constant inlet energy level but is increasing. 


Perhaps the major discrepancy between the experimental 


test results and the computer analytic model is evidenced 
in the axial pressure profiles of Figures 36 and 37. It is 
clearly evident that at a certain point the nozzle passage 
reverts from a nozzle to a diffuser. This transition point 
occurs slightly downstream of the throat and the pressure 
starts to increase. 

Unfortunately the computer analytic model requires a 
pressure profile as an input. Furthermore the pressure pro- 
file must be continuously decreasing. Thus the profile as 
obtained from the experiment are not directly useable. The 
situation is examined with the aid of Figure 38. Curve A is 
a typical axial pressure distribution as obtained from the 
experiment. Curve B is a pressure distribution used by 
Elliott in the application of his computer program. Pressure 
profiles C and D were a bitrarily defined and the exit plane 
velocity for each was calculated. Velocity variation in the 
range of 1% is evident. It appears that the final exit 
velocity is relatively insensitive to the actual pressure 
crofile in the nozzle. 

In light of the preceeding difficulty cf matching the 


Y 


experimental pressure profile te the computer mcdel and the 
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apparent flexibility of the type of pressure profile it was 


decided to employ a form of profile D of Figure 38. Thus 
instead of using the actual experimental pressure data as 
an input in the computer model, a profile resembling curve 
D was developed for each inlet pressure case. 

It is apparent that within certain rather wide limits 
the Elliott computer model yields results that correspond 
within 10% to results obtained from the experiment. The 
general trends have been confirmed and their behavior has 
revealed nothing unexpected. The conclusion of Elliott 
C[Ref. 1] has been largely substantiated. "It is very 


difficult to design a bad or a good dual-phase nozzle." 
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Figure 26. Mass Ratio vs. Exit Velocity at 
Pressure = 29+1.S5psi Experimental 
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Figure 27. Mass Ratio vs. Exit Velocity at 
Pressure = 35+1.5psi Theoretical 
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Figure 78. Mass Ratio vs. Exit Velocity at 
Pressure = 29+1.S5psi Theoretical 
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Figure 30. Mass Ratio vs. Exit Velocity at 
Pressure = 29+1.5psi Area = .842375 
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Figure 32. Mass Ratio vs. Exit 
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VIII. CONCLUSION 


A set of experiments were conducted with an air-water 
mixture flowing at a mass ratio from 2 to 13. The experi- 
ments correlated well with a dual-phase two-component 
computer program. It appears that the program will permit 
prediction of the exit velocities to an accuracy of 10% 
for two-phase nozzles. At small nozzle area ratios, accu- 
racies between theoretical and experimental are better; 
approximately 5%. All predicted velocities are higher 
than measured experimental values due to estimation of the 
drop size and the initial kinetic energy of the liquid at 


: 


njection point. 


& 


b+ 


It would be desirable to develop a drop size subroutine 
to better estimate varying drop size. The experimental 
system can be improved by velocity measurement devices at 


the inlet. Input of the gas and liquid velocities are 


necessary in the dual-phase two~component program and thus 
these measurement devices are vital for better accuracy. 

Due to the insensitivity of the pressure distribution, 
it appears that any reasonable approximation to the pressure 
profile can be employed in the two-phase two-component flow 
program. For a given nozzle exit area and inlet pressure, 

a pressure distribution curve can be obtained threugh the 
experimental system. This distribution can therefore be 
used to obtain good results. 
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APPENDIX I: CALIBRATION PLOT OF PRESSURE TRANSDUCERS | 
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APPENDIX J: CALIBRATION PLOT FOR FORCE-BLOCK 


Calibration Plot For Force-Black 
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PRINT “INPUT AOTUMETER READING” 


3 

4 

5 PRINT Rr 
7 DIM A(14),B(14),C(14) Dist(14), D(14) 

8 DIM Press(14) 

10 rt X(18,100).Sum( 18) ,Mean(18), Sum2(18) 


N+10 
40 FOR I=1 TON 


S50 OUTPUT 709;"AI1VT1" 
60 ENTER 709:X(1, T) 
70 OUTPUT 709; ral2vt1* 
30 ENTER 709:X(2,1I 
30 OUTPUT 709: sorauTte 
100 ENTER 709:X(3.1) 
110 OUTPUT 703; “ALGVTI" 
120 ENTER 709:X(4,1) 
130 OUTPUT 709; “atsytt" 
140 ENTER 709: X(5,1 
150 QUTPUT 709; “alevTi" 
160 ENTER 709:X(6,1) 
170 QUTPUT 709; ratzvri* 
180 ENTER 709; X(7,1 
190 OUTPUT 709; ‘aigurt* 
2090 ENTER 709; x(8, 
210 OUTPUT OF et AtovTa* 
220 ENTER 709:X(9,1) 
230 OUTPUT 769: “ALTOVTI® 
240 ENTER 709:X(10.1) 
250 OUTPUT 709;“AI1TIVTI" 
260 ENTER 709; XC14, 1) 
261 OUTPUT 709; “*AIT2VTI" 
262 ENTER 709; X(12,1) 
263 OUTPUT 709;"AI13VT1)° 
264 ENTER 709:X(13,1) 
265 OUTPUT 709: "AL14VT1)° 
266 ENTER 709; X(14,T) 
267 OUTPUT 709: "AIISVT1)" 
268 ENTER 709;X(15,1) 
270 NEXT 
440 FOR J=1 TO 16 
450 Sum( J) =0 
460 FOR I*1 TON 
470 Sum(J)*Sum( J) +X J. 1) 
480 NEXT I 
490 Mean(J)=Sum(J)/N 
510 NEXT J 
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